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ABSTRACT 



Aims. We present the results of a spectroscopic search for close binaries among horizontal branch (HB) stars in NGC6752. 
Methods. We used the ESO VLT-FORS2 instrument to obtain medium resolution (R=4100) spectra of 51 hot HB stars with 8 000 K < Teff < 
32 000 K during four consecutive nights. Eighteen of our targets are extreme horizontal branch (EHB) stars with Teff > 22 000 K. Radial 
velocity variations were measured with cross-correlation techniques and we carefully evaluated the statistical and systematic errors associated 
with them. 

Results. No close binary system has been detected among our 51 targets. The data corrected for instrumental effects indicate that the radial 
velocity variations are always below a 15 km s"' (3cr level). From a statistical analysis of our results, we conclude that (at 95% confidence 
level) the fraction of binaries with a ~ 0.5 Mq companion among EHB stars in NGC6752 is smaller than 20%. 

Conclusions. This empirical evidence sharply contrasts with what has been found for hot subdwarfs in the field, and opens new questions about 
the formation of EHB stars in globular clusters (and possibly in the field as well). 

Key words, stars: horizontal branch - binaries: close - binaries, spectroscopic - globular cluster: individual: NGC6752 
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1. Introduction 

Although stellar evolution theory has successfully identified 
horizontal branch (HB) stars as post-helium flash stars of 
low ( typically 0.7 - 0.9 Mp in^ globular cluster s) initial 
mass dHovle & Schwarzschild"l955; "Faulkner"l966'), we still 
lack a comprehensive understanding of thei r nature. The large 
morphological differences (see for example IPiotto et alJl2002l) 
among HBs of Galactic globular clusters (GCs) are surely 
one of the most puzzling problems of stellar evolution. It is 
well known that cluster metall icity affects the HB morphology 
dSandage & W allerstein"l960'). but there are still many mor- 
phological differences among the HBs of GCs with the sam e 
metallicity dSandage & WildevI 1 19671 Ivan den Berg hl Il967t) . 
suggesting the presence of one or more second parameters. 
Many "second parameter" candidates have been proposed, but 
none has provided an overall explanation for all the available 
observations, and it is likely that the so-called "second parame- 
ter problem" is the cons equence of a combination of parameters 
jFiisiPecci etal.ll993l) . 

In more recent years, more and more GCs have been found 
to show an HB blue tail extending toward increasing effec- 



* Based on observations with the ESO Very Large Telescope at 
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live temperatures and fainter visual magnitudes ("Sosin et al.| 
'l997; Rich et al. 1997; Ferraro et al. 1998; Piotto et al. 1999). 
This extreme horizontal branch (EHB) population with Teff > 
20 000 K is a challenging test of the theoretical models of 
the late evolutionary stages of low-mass stars, and their ori- 
gin is still not understood (see below). The EHBs might also 
have important implications for extragalactic astronomy, since 
EHB stars have been identified as possibly being responsible 
for the UV upturn in elliptical galaxies. This UV upturn cor- 
responds to the increase in flux with decreasing wav elengt h 
below 2500 A discovered in the bulge of M31 (Cod3 ll969t) . 
and then found in almost all bright elliptical galaxies and 
many spiral bulges. It has been proposed that the EHB stars 
are the so urce o f this in crease in UV flux (see f or example 
ICreggio & Renzinilfl 99d 11995 iBrown et~Zil200nl) . This idea 
has been reinforced by the discovery of EHB stars in metal- 
rich GCs like NGC6388 and NGC6441 (Rich et al. 1997) and 
in the Galactic bulge dBusso et alJ2005ll) . As EHB stars evolve 
with time, the upturn is expected to change during the evolution 
of the galaxies, and to be almost absent for younger ellipticals 
dGreggio & Renziniill99d iTantalo et alJll996l) . A knowledge 
of the formation mechanisms and evolution of EHB stars in 
different environments is therefore urgently needed in order to 
model the behavior of the UV upturn with age (redshift). The 
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UV upturn could then also potentially be used as an age dis- 
criminator for elliptical galaxies. 

While it is agreed that EHB stars are He-burning stars 
that have suffered heavy mass loss during their ev olution 
lllben & Roodll970HFaulkneill972HD'Cruz et al.ll996l) . keep- 
ing only a thin external envelope with a mass around 0.02 M0 
or lower, their actual formation mechanism remains un- 
clear. Some mechanism able to enhance the mass loss must 
be r esponsible for their formation. As first proposed by 
IWils on & Bowen ( 1984), many authors explored mass-loss 
mechanisms of the HB itself. lYong et aO (l2000l) successfully 
reproduce the morphology of blue HBs in metal-rich GCs with 
a constant mass-loss rate during HB evolution, b ut provide 
no explanation for the required rate. IVink & Cassisi 120021) 
show that enhanced radiation-driven winds cannot account for 
the high mass-loss rates required, concluding that the most 
plausible mechanisms are unable to form EHB stars by mass 
loss on the HB. On the other hand, many authors propose 
heavy mass-loss rates during the previou s red giant branch 
(RGB) phase (see for example ISoker et al .' 2001). IbenI ( Il99d) 
propose an alternative scenario, in which the EHB stars are 
formed from mergers o f helium white-dwarf binary systems. 
iFusi Pecci et indicate that the GC density can en- 

hance the presence of EHB stars, hinting that dynamical in- 
teractions may play an important role in their formation. 

There is still a lack of models that can explain such heavy 
mass loss in the evolution of a single low-mass star. The bina- 
ritv of EHB stars, as pr oposed by many authors ( Mengel et al. 
ll976tlHeber et alJ2002; '). might provide an explanation for their 
formation, since the dynamical interaction with a close com- 
pact companion can enhance the mas s loss through a number 
of different binary-evolution channels dHan etal.l2002i see be- 
low), particularly during the RGB phase. The pres ence of a 
binary population in GCs is now well established llHut et alJ 
Il992: Bailvn' l993h . and binaries can get closer and cl oser as a 
consequence of the dynamical interactions inside GCs jHeggid 
Il975h . The idea that EHB stars are components of close binary 
systems has been strengthened by observations of field EHB 
stars, also known as subdwarf B-type stars (sdBs). Binaries 
have b een found to be very common among them. Maxted et al. 
( I2OOII) conclude f rom their observations that 6 9 + 9% of sdB 
stars are binaries. iMorales-Rueda et"aD ( l2003l) recently mea- 
sured the orbital periods P and the semiamplitudes of radial 
velocity variation K of 22 new binary sdBs, increasing the 
number of sdBs for which these quantities are known to 38. 
Thirty of them have periods below 3 days (22 with P < 1 day), 
showing that they are almost all very close systems, with K 
usually exceeding 50 km s"' (only 7 exceptions) and easily 
greater than 100 km s ~'. A more recent sear ch for binaries 
among field sdB stars (Na oiwotzki et alJ l2004') led to a much 
lower (42%) close bin ary fra ction than expected from previ- 
ous results. iNapiwotzk i et alJ (|2io04) point out that their sam- 
ple contained a much higher percentage of faint, i.e. distant, 
stars than the one of Maxted et al.. (.200 U . The Napiwotzki et 
al. sample may therefore extend to greater distances from the 
Galactic plane and thus be contaminated by thick-disk or halo 
members. This could imply that a relation of binary frequency 
with metalUcity and/or age is present. 



iHan etalJ tOO^ find from binary population synthesis 
technique s that 7 6-89% of the sdBs should be close binaries. 
iHan et all ll2002 ) analyzed in detail the main binary evolution 
channels that can lead to the sdB formation. Comparing their 
models with the available empirical data, they find that a very 
efficient common-envelope (CE) channel can fit the observed 
distribution of close-binary sdB periods, though it is not pos- 
sible to fi t the observed cum ulative luminosity functions, as 
shown by Lisker etaP ( l2005l) . showing that there are still some 
problems with the binary scenario. In the CE channel, the pro- 
genitor of the sdB star is a giant that fills its Roche lobe near the 
tip of the red giant branch and experiences mass transfer dur- 
ing which the core of the giant and the companion (a normal 
main-sequence star or a white dwarf) spiral towards each other 
inside a common envelope formed out of the giant envelope. 
Once enough orbital energy has been released, the envelope is 
ejected and the system becomes a short-period binary. 

There are alternative mechanisms suggested for the for- 
mation of t he EHB stars. Recent observational results on 
O) Centaur i ('Bedinetal. 20o3 IPiotto et alJ boOSi). M3-M 1 3 
(ICal oi & D'Antona .2005 ). and NGC2808 jP' Antona eTall 
2005) have given new impulse to the idea that EHB stars 
may result from a second generation of stars enriched in he- 
lium by pollution from intermediate-mass AGB star ejecta 
(D'Antona et al. 2002; Lee et al. 2005). 

In the attempt to shed some light on the origin of EHB stars 
in globular clusters, we want to test in this paper the possibil- 
ity that they are stars that have experienced greatly enhanced 
mass loss during the evolution in a close binary system. We 
started our spectroscopic search for close binaries among the 
EHB stars of NGC6752. This GC is an ideal target as it is dy- 
namically evolved, shows an extended and well-populated blue 
HB (as shown in Fig. a large population of X-ray sources 
in the core that could be cataclysmic variables (Poolev e t alJ 
l2002h . and it is supposed to have a large main-sequence binary 
population jRubenstein &Bailvnll997.) . 

In a forthcoming paper, we will present and discuss the 
main atmospheric parameters (Teff, iogg, log mass) derived 
from a complementary set of spectra obtained during the same 
observational run. In the following, we concentrate on the ra- 
dial velocity and radial velocity variation measurements, and 
on their implication for the close binary scenario in the forma- 
tion of EHB stars. 

2. Observation and data reduction 

The spectra were acquired during four nights of observation, 
from June 11 to June 14, 2002, at the VLT-UT4 telescope 
equipped with the FORS2 spectrograph in MXU mode. Fifty- 
one HB stars were selected from th e photometric catalog of 
NGC6752 by Momanv et alJ 1I2OO2I1 . The stars are well dis- 
tributed along the HB from the cool edge (Tetr ~ 8 000 K) to 
the EHB (Teff ~ 30000 K), and were divided into 3 fields for 
multi-object spectroscopy. The position of our targets on the 
HB are indicated in Fig.[2and their radial distribution is shown 
in Fig. 2. 
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Table 1. List of the observed stars wi th their physical parameters. IDs, coordinates (RA and DEC), magnitudes V and color 
(U - V) are from lMomanv et al] (l2002l) . The absolute heliocentric radial velocities (Vrad) were measured in the present work, the 
effective temperatures were obtained from the U - V color of the targets, with the relation (U - V) vs. Tgff shown in Fig^and 
described in the text. The estimated error for temperatures is o-j- 5% 



field slit 


ID RA (J2000) 


DEC (J2000) 


V 




Teff (U-V) 




hh:mm:ss 






km s"' 


K 



A 


1 


14770 


19 


11:32.563 


-59:59:37.63 


17.247 


-35 ± 


10 


30300 


-1.192 


A 


2 


11634 


19 


11:16.084 


-60:00:27.33 


14.431 


-43 ± 


10 


10000 


-0.174 


A 


3 


14944 


19 


11:22.485 


-59:59:35.51 


15.401 


-44 ± 


10 


15600 


-0.671 


A 


4 


15026 


19 


11:12.319 


-59:59:34.51 


14.018 


-36 


± 9 


9100 


-0.042 


A 


5 


16551 


19 


11:14.253 


-59:59:13.90 


15.420 


-44 ± 


10 


15500 


-0.663 


A 


6 


15395 


19 


11:05.940 


-59:59:29.51 


17.337 


-40 


± 8 


29000 


-1.153 


A 


7 


20919 


19 


11:10.496 


-59:58:15.25 


13.844 


-34 


+ 7 


8100 


+0.110 


A 


8 


18782 


19 


11:16.428 


-59:57:45.90 


14.601 


-35 ± 


10 


13100 


-0.498 


A 


9 


17941 


19 


11:04.827 


-59:58:02.02 


15.928 


-39 


± 8 


16700 


-0.674 


A 


10 


20302 


19 


11:16.912 


-59:57:13.20 


16.594 


-36 


± 7 


20800 


-0.901 


A 


11 


26756 


19 


11:25.759 


-59:56:17.07 


14.599 


-21 


± 9 


10800 


-0.265 


A 


12 


27181 


19 


11:24.177 


-59:56:04.10 


15.073 


-34 ± 


10 


15400 


-0.661 


A 


13 


24849 


19 


10:59.399 


-59:57:05.99 


14.578 


-49 


± 9 


10800 


-0.273 


A 


14 


27604 


19 


11:17.021 


-59:55:50.56 


15.937 


-36 


± 9 


19200 


-0.841 


A 


15 


28231 


19 


11:18.560 


-59:55:27.51 


17.367 


-29 


± 6 


26900 


-1.091 


A 


16 


26760 


19 


10:59.788 


-59:56:17.51 


15.515 


-37 


± 7 


16400 


-0.716 


A 


17 


28554 


19 


11:07.939 


-59:55:14.16 


17.289 


-19 


± 8 


27600 


-1.114 


A 


18 


28693 


19 


11:08.108 


-59:55:08.53 


17.306 


-26 ± 


10 


29600 


-1.172 


A 


19 


28947 


19 


11:04.724 


-59:54:56.61 


16.832 


-36 


± 8 


23200 


-0.979 


B 


1 


4964 


19 


10:49.893 


-60:04:08.04 


14.552 


-47 


± 8 


10700 


-0.258 


B 


2 


49317 


19 


11:15.879 


-60:06:00.17 


13.889 


-35 


± 8 


8600 


+0.028 


B 


3 


5455 


19 


10:52.734 


-60:03:36.84 


17.371 


-34 


± 6 


28700 


-1.147 


B 


4 


5487 


19 


10:57.554 


-60:03:35.02 


16.766 


-43 ± 


10 


21800 


-0.935 


B 


5 


5134 


19 


11:09.624 


-60:03:56.23 


15.615 


-40 


± 8 


16200 


-0.706 


B 


6 


4672 


19 


11:21.548 


-60:04:31.46 


17.181 


-33 


± 8 


27900 


-1.123 


B 


7 


5201 


19 


11:26.748 


-60:03:51.85 


17.615 


-37 


+ 7 


29900 


-1.181 


B 


8 


5865 


19 


11:21.577 


-60:03:14.90 


17.195 


-44 


±9 


25100 


-1.037 


B 


9 


7843 


19 


11:14.678 


-60:021:55.92 


15.317 


-41 


±8 


16000 


-0.693 


B 


10 


6284 


19 


11:29.780 


-60:02:54.90 


17.155 


-36 


±6 


28600 


-1.142 


B 


11 


10257 


19 


11:07.698 


-60:00:53.00 


14.075 


-37 ± 


10 


9300 


-0.075 


B 


12 


10625 


19 


11:21.006 


-60:00:45.37 


17.701 


-24 


±9 


30700 


-1.201 


B 


13 


8672 


19 


11:32.600 


-60:01:30.83 


17.790 


-33 


±9 


33500 


-1.266 


B 


14 


10711 


19 


11:28.425 


-60:00:43.55 


17.418 


-31 


±7 


30200 


-1.189 


c 


1 


11609 


19 


10:36.021 


-60:00:28.10 


15.311 


-37 


±9 


14600 


-0.613 


C 


2 


14664 


19 


10:41.367 


-59:59:39.90 


13.972 


-7 


±9 


8300 


+0.086 


C 


3 


14727 


19 


10:39.624 


-59:59:39.13 


14.253 


-33 


±9 


9600 


-0.119 


C 


4 


35186 


19 


10:13.293 


-60:00:03.25 


14.443 


-34 


±9 


10600 


-0.239 


C 


5 


35662 


19 


10:22.119 


-59:59:30.30 


15.010 


-40 


±9 


14000 


-0.566 


C 


6 


35499 


19 


10:10.920 


-59:59:41.57 


14.895 


-35 


±9 


13400 


-0.518 


C 


7 


36242 


19 


10:23.020 


-59:58:47.17 


15.024 


-31 


±7 


14100 


-0.570 


C 


8 


36480 


19 


10:23.013 


-59:58:30.37 


16.899 


-35 


±8 


24800 


-1.028 


C 


9 


36502 


19 


10:18.896 


-59:58:28.74 


14.863 


-35 


±9 


13200 


-0.505 


C 


10 


36830 


19 


10:02.476 


-59:58:03.50 


17.388 


-33 + 8 


30700 


-1.202 


C 


11 


38095 


19 


10:26.616 


-59:56:23.29 


15.358 


-35 


±9 


16100 


-0.700 


C 


12 


38087 


19 


10:18.120 


-59:56:23.89 


17.190 


-37 


±8 


30500 


-1.197 


C 


13 


32470 


19 


10:22.172 


-59:55:54.36 


14.681 


-27 


±8 


10800 


-0.271 


C 


14 


28695 


19 


10:35.134 


-59:55:08.70 


14.329 






9900 


-0.155 


C 


15 


38504 


19 


10:15.308 


-59:55:41.62 


14.823 


-39 


±9 


12900 


-0.155 


C 


16 


39008 


19 


10:25.034 


-59:54:41.75 


17.178 


-15 ± 


11 


31400 


-1.219 


c 


17 


38889 


19 


10:09.836 


-59:54:56.10 


14.897 


-48 


±9 


13400 


-0.520 


c 


18 


38963 


19 


10:03.450 


-59:54:46.56 


14.264 


-17 ± 


10 


9600 


-0.110 
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Fig. 1. Positions of target stars (solid dots) in the color- 
magnitude diagram of NGC6752. Data from Momanv et al. 




Fig. 2. Radial distribution of the observed stars. The dark 
shaded area indicates hot stars (Teff > 20000 K). The half-mass 
radius from Harris (1996) is also indicated. 



The list of targets, their coordinates, V magnitudes, and the 
relevant parameters extracted from the color-magnitude dia- 
gram and our spectra are given in Table The approximate 
effective temperatures were obtained from the U -V colors by 
iMomanv et alJ ll2002h . by applying the color-temperature rela- 
tio n derived from the te mperatures spectroscopically measured 
bv lMoehler et alJ 12000) on a sample of HB stars in the same 
cluster as in Fig. |3] Combining the dispersion of the points 
around the fitted line in Fig. |3j and the photometric errors in 
U - V we estimated the error in effective temperature to be 
approximately 5%. 

During each night, up to two pairs of 1 800s exposures were 
secured in each field with grism 1400V-I-18, as shown in Table 
|2] The slit width was O'.'S and the resulting resolution 1.2 A. 
The exposures were always collected in pairs and subsequently 
summed (see ^3.1> except on the 3rd night for field A, when 
only one 1800s exposure was acquired. Just before each pair 
of exposures, a slit image (without grism) was taken, and these 
frames were used in the correction of instrumental effects that 
affected the data ( ^3.2\ . The bias, flat, and lamp images were 
acquired before and at the end of each observing night. The 
observed spectral range was a: 1300 A wide but, due to the 
different positions of the slits within the masks, the central 
wavelengths of the spectra varied. The bluer spectra reached 
X! 4200 A in the blue edge encompassing the Hy line and the 
447 1 A Hel doublet, while others extended on the red side up 
to X 6000 A. The line was always inside the observed spec- 
tral range, except for star 28695. This cool target was excluded 
from our analysis due to the lack of strong lines in the spectral 
range. 



Data reduction from the multi-spectrum frames to the one- 
dimensional calibrated spectra was performed with standard 
MIDAS' procedures. All slitlets were extracted from the full 
frames (bias, wavelength calibration, flat field, science spectra) 
and reduced independently. The wavelength calibration (wlc) 
was performed with the HeNeHgCd lamp images, fitting a 3'^'' 
order polynomial to the dispersion relation. We re-binned the 
2D frames to constant wavelength steps of 0.25 A/pix. The ex- 
traction region for target and sky varied strongly from star to 
star because different parameters affected the selection, mainly 
the crowding level and the object magnitude. We tried to ex- 
tract all the spectra with an optimal extraction routine jHornd 
'1986'), but for some stars (mainly the brightest ones) the proce- 
dure gave bad results, so in these cases we opted for a simple 
average over the extraction region. 

No flux calibration was applied to the spectra, because they 
were continuum-normalized for the measurement of radial- 
velocity variations. 

3. Measures 

3.1. RV-variation measurement 

Radial velocities (RVs) were measured wit h the cross- 
correlation (CC) technique jTonrv & DaviJ 19791) . using the_^- 



' ESO-MIDAS is the acronym for the European Southern 
Observatory Munich Image Data Analysis System that is de- 
veloped and maintained by the European Southern Observatory 
l_http://www.eso.org/projects/esomidas/j 
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Table 2. List of 1800s exposures acquired each night. The UT 
of the start of the exposures (hour and minutes) is indicated, 
and also the averaged seeing (variations between the exposures 
in the same night are within O'.'l). 



field 


night 




12 


13 


14 


15 






2:44 






A 


8:43 


3:15 


9:34 


7:47 




9:14 


3:56 
4:27 


8:18 






5:55 








B 


6:27 


5:08 


8:33 


8:58 




6:59 


5:39 


8:59 


9:29 




7:30 












6:39 






C 




8:04 


7:05 


6:36 






8:34 


7:33 
7:59 


7:06 


Average seeing 


a: 9 


1'.'2 


1'.'3 


1'.'4 



cor IRAF^ task. Before co-adding each pair of spectra we mea- 
sured the RV variation between them, in order to verify that no 
significant variation had occurred and that no information was 
lost due to averaging. The summed spectra will be referred to 
with the field and the night of observation (from 12 to 15), with 
an additional letter (a or b) when two pairs were acquired in 
the same night and field. The A14 spectra come from a single 
1 800s exposure, since only one frame was acquired that night, 
and in our analysis we took their lower S/N into account. 

For each target star all the spectra were cross-correlated, 
thus performing 10 CCs for each star in fields A and B and 6 in 
field C, covering different temporal intervals from one hour to 
3.1 days. In each measurement the first spectrum (in temporal 
order) was assumed as template (the choice of which of the 
two spectra is used as template affects only the sign of the RV 
variation). 

For each RV measurement the cross-correlation function 
(CCF) was computed and the positi on of its peak determined 
with a Gaussian fit (see for examplel Recio-Blanco et alJl2004t 
IPubath et al. 1990, for a description of the procedure). In 
the measurement concerning Balmer lines, the CCFs have 
Gaussian shapes only in the peak but not in the wings, and 
sometimes they are asymmetric; whereas the CCFs obtained 
from weak lines are approximated well by a Gaussian, but they 
are often distorted due to low S/N. Therefore the fit was re- 
stricted to the central portion around the peak, but different fits 
were tried, varying the width of the CCF region fitted, in order 
to find the best estimate for the RV variations. 

Our analysis focused on the Hy? line, and 432 measurements 
were performed in the 4 830 - 4 890 A spectral range (H^ line 
with full wings). We also compared these results with measure- 
ments from other parts of the spectrum, both to clarify the influ- 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 




-0.5 

U - V 



Fig. 3. Plot used to derive the approximate eff'ective tem- 
peratures of our target s. The poin ts are the eflFective tem- 
peratures measured by iMoehler'et alJ ( l2000l) using models 
with solar metallicity, plotted against the U - V color in 
the Momanv et al. (2002) photometry. Seven stars listed in 
(Moehler et alJ ll2000l) have been e xcluded since their identifi- 
cation in the lMomanv et al] ( l2002h photometrv is dubious. The 
line indicates the Teff vs. color relation obtained by a fit of the 
data points, used to estimate the temperatures of the stars ob- 
served in our work. The point at f/ - V ^-0.8, Teff ~ 29000 K 
were excluded from the fit. 



ence of noise on the measured RV variations and to overcome 
the intrinsic weakness of the CC tech nique for blended line s 
in double spectra (see for example IZucker & Mazehl|l994h . 
Therefore, although our survey was intended to be based on 
analysis of the H^, all the measurements were repeated, cross- 
correlating the entire spectra (with and without H/j) and Hy 
(alone and with H^) when it fell inside the spectral range, for 
a total number of 1 532 CCs. The presence of many metal- 
lic lines in the spectra of stars with 11 500 < Teff < 18 000 
K (see Fig. Bl, due to radiative levitation of heavy elements 
dCilaspev et aDll989HBehJl2003l) . gave good CCFs even with- 
out the H lines. On the other hand, the low S/N and the lack 
of useful lines in the spectra of hot stars (Teff > 18,000K) 
usually prevented cross-correlating the entire spectra without 
H^, and then only spectral sections with the strongest helium 
lines were cross-correlated in place of these measurements, 
with quite uncertain results (see §14. It . Fourier filters of vari- 
ous shapes jBrault & Whitelll97ll) were applied to almost all 
the noisy spectra, obtaining better CCFs but unchanged results 
for the most part. 

The [OI] 5577 A sky line was used as zero-point in order 
to correct possible spectral shifts due to differences between 
lamp and star spectra. The sky line fell outside the spectral 
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Fig. 4. A sample of the collected spectra, normalized and ordered by increasing Teff . The increasing noise with temperature is 
clearly visible, as well as the large quantity of metallic lines for stars with 1 1 500 < Tefi < 18 000. 



range for 13 stars, preventing us from correcting the spectra 
before the measurements. We found that the differences in the 
sky line position in the spectra between each pair of frames is 
a linear function of Y-position (perpendicular to the dispersion 
direction) of the slits in the masks (Fig.|5}. From the measured 
sky line positions, we then calculated these differences for each 
pair of spectra, and applied them as corrections to the RV vari- 
ations measured with the CCs. For the 13 stars missing the [OI] 
line, we calculated the corrections from their Y-coordinates by 
means of the least-square fit obtained by using the spectra with 
the [OI] fine. 

3.2. Correction of the RV measures 

After the skyline-based correction mentioned above, we were 
forced to correct the measured RV variations again, since they 
showed a clear correlation with the displacement (along the 
dispersion direction) of the stars with respect to the center of 
the slit. The eff'ect resulted in an evident systematic error, up 
to 10-12 km s ^ of the same order of magnitude for all the 
stars within the same pair of frames (usually with a dependence 
on Y-position on CCD, probably due to rotation of the mask 
and/or the field). In the slit frames we measured the position 
of the stars relative to the center of the slits with a Gaussian fit 
of the stellar profile parallel to the dispersion direction. Then, 
we obtained the differences in these positions between pairs 
of frames, and translated them from pixels to km s ' with the 
instrumental relation 38.2 km s ' pixel in order to evaluate 
the effect of this shift on the RV measurements. Actually, it is 
not the real stellar profile that has been fitted, but the star pro- 
file convolved with the narrow slit, that generates the spectrum 
on the CCD. In fact, for narrow slits, the relation between the 



movement of the star behind them and the induced RV vari- 
ation should be flatter than the dispersion relation, but a shift 
of the center of the profile that reaches the grism through the 
slit is expected to induce an RV variation given by the disper- 
sion relation, with good approximation. This is confirmed by 
the plot in Fig.|6l where we compared the two quantities (shift 
of the profile and RV variation) for each star and each pair of 
frames: the points show a good agreement with the expected 
relation indicated by the straight line. Although the plot con- 
firms that it is a good solution for relating the two quantities, it 
also possibly indicates that it is just a first-order approximation, 
because a deviation from the straight line for higher displace- 
ments can be seen. The real solution could be an S-shaped func- 
tion. Nevertheless we adopted a linear approximation, that stiU 
gave good enough corrections for our purposes and our overall 
errors, in order to correct the RV variations with quantities that 
depend only on the measured displacements and not on the RV 
variations themselves (as would be the case if, for example, we 
fit the points in Fig.|6lto obtain a higher-order relation). 

We then compared the displacements of the profiles and the 
RV variations for each star, as a function of the Y coordinate 
on CCD. A typical example is shown in Fig. (upper panel). 
A dependence on slits' Y-position can also be seen, but with a 
certain scatter due to the random errors that affected the stellar 
fitting procedure in the narrow slits. We preferred to derive the 
final corrections by applying the value of the least-square fit of 
these plots to each star instead of the real shifts measured star- 
by-star, in order to avoid introducing additional noise into the 
data. 

For frames B12b and A 13a, the slit image was not ac- 
quired immediately before the exposures, but one hour before 
for frame B12b and one hour after for A13a. In both cases, 
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Y coordinate on CCD (pixels) 

Fig. 5. An example of the corrections to re-position the 5577 A sky Hne to its laboratory wavelength. The differences (in km s"') 
in the position of the sky line in the spectra A15 and A12 are plotted as a function of the Y-positions of the slits on the CCD (in 
pixels). The corrections for the few stars without the sky line in the observed spectral range hwere obtained from their Y-position 
by the least-square solution (straight line in the figure). The error bar in the plot is the Icr estimated for this corrections (see Sect. 
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Fig. 6. Displacements of stars with respect to the center of the 
slit plotted against the RV variations measured from CCs. The 
straight line indicates the instrumental relation 1 pixel = 38.2 
km s"' given by the dispersion of the spectra. 



the corrections do not agree with the RV variations observed, 
but a constant shift between them is clear (Fig. middle 
panel) and has of the same amount for both frames. This in- 
dicates a time-dependent instrumental movement (probably a 
shift of the mask in its housing) that shifted the spectra slightly 
7 km s"' = 0.2 pixels). Either way, the slit image for frames 
B12a and B12b was the same (collected just before the first, 
B12a), and they were taken in sequence, without any telescope- 
pointing change, and calibrated with the same lamp image. For 
all these reasons, it is reasonable to assume that the observed 
shift of the sky Une on the CCD between the frames B12a and 



B12b was only due to this time-dependent instrumental move- 
ment. Therefore we think that we have been able to measure 
and correct this effect properly. 

The corrections derived for frame B13 do not overlap the 
RV observed variations (Fig.0 lower panel). Since it is the only 
frame out of 14 for which this happens, its slit frame seems not 
to be reliable for some reason. Probably the mask moved in- 
side its frame between the slit image acquisition and the spec- 
tra exposure. We estimated a mask shift of 0.13-0.34 pixels 
with a rotation 0.00928 degrees. This is not unlikely and simi- 
lar effects have been noted by other observers before. Since we 
could not both find an explanation and prove it, we preferred to 
exclude this frame from further analysis. The uncorrected RV 
variations observed are still useful and give some information 
(see although only in a qualitative way. 

3.3. Absolute RV measures 

We measured absolute RVs in order to check the cluster mem- 
bership of the observed stars by means of CCs with the tem- 
plate star HD1881 12, a binary sdB star with known ephemeris 
( Heber et al.,,2QQ3). The results are shown in Col. 7 of Tabled 
The final absolute RVs have been corrected for sky line posi- 
tion, as already described in ^3.21 All the stars show an abso- 
lute RV in a greement with that of the cluster (-27.9 km s~\ 
iHarrisI 19961) within 2cr, and can thus be considered RV cluster 
members. 



4. Error analysis 

We performed an accurate error analysis in order to estimate 
the significance of the measured RV variations. We analyzed 
all the main error sources, and estimated a standard error cr for 
each of them. The final errors associated to the results presented 
in Table|4]were obtained as their quadratic sum. 
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Fig. 8. Shift of lamp lines compared to their laboratory wavelength, in km s ' . Each circle is proportional to the frequency of the 
datum. The dashed lines indicate the Xcr dispersion. 
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Fig. 7. Comparison between the measured displacement of 
stars inside the slit (solid points), translated from pixels to km 
s~\ and the RV observed variations (open points). The lines are 
the least-square solutions for the two sets of data. Upper panel: 
a typical plot (frames A15 and A 12), in which the two sets of 
data are very similar Central Panel: same plot involving frame 
B12b. A shift between the two lines is evident. Bottom panel: 
same plot involving frame B13. The two sets of data do not 
agree. 



4.1. Cross-correlation errors 

The error in the CC procedure is obtained by the CC theory as 



mi + r) 



where 



V2cra 



and N is the number of bins, B the highest wavelength in which 
the CCF Fourier transform has an appreciable value, h is the 
height of the CCF peak, and the rms of the antisymmetric 



part of the CCF (see lTonrv & Davisll979t for details). This er- 
ror is computed directly by thtfxcor IRAF task during each CC 
procedure. In the measurements with Balmer lines it is usually 
between 0.5 km s"' (for high S/N spectra) and 2 km s"' (for 
asymmetric CCFs), but occasionally CCFs with bad profiles 
lead to a CC error up to 5 km s The measurements without 
Balmer lines are affected by much higher errors for hot stars 
(usually 3-5 km s"'), whereas for cooler stars the accuracy re- 
mained quite unchanged due to the numerous metallic lines in 
their spectra. 

4.2. Wavelength calibration errors 

The precision and reliability of the wlc procedure, which is es- 
sential for our analysis, was tested in each of the 186 lamp im- 
ages used. These images were calibrated with the coefficients 
obtained in the wlc procedure, and the position of 9 bright lamp 
lines were measured with a Gaussian fit and compared to their 
laboratory wavelength. This analysis pointed out the level of 
inaccuracy in the wlc procedures that could affect the RV mea- 
sures. The obtained data are plotted in Fig.|8] 

No systematic error was found. Only the lines 4471 and 
5875 A seem calibrated redder than the theoretical wavelength, 
but every night the mean shift of the former was not greater 
than the wlc error, and the latter line lies in the extreme red part 
of the spectra, never used in the CCs, and between two well- 
calibrated lamp lines. The dispersion of the data is in the range 
1.3-1.7 km s"', slightly but not significantly varying among the 
nights and the three fields. They were considered as an estimate 
of the wlc error oSmXc, which has been counted twice (in quadra- 
ture) in the final error estimate because in each CC two spectra 
are involved. 

4.3. Extraction and fit errors 

Two additional sources of error are present in the measurement 
procedure. The first one is related to the extraction of spectra 
from the 2D images, as the choice of the aperture width can 
change the profile of spectral lines (mostly in noisy spectra or in 
the presence of cosmic spikes) slightly, and the CC procedure 
is sensitive enough to reveal the differences. The second, more 
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Table 3. Standard deviations (km s"') of the differences be- 
tween the RV measurements obtained with different fits and 
spectral extractions In field B the bluest stars have been dis- 
tinguished from the others, and "single spectrum" indicates the 
measurements involving the not-summed spectra A14. 
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important, is related to the choice of the number of CCF points 
used to compute its Gaussian fit. For distorted CCFs, this eiTor 
is much greater than crcc, since the final result could depend 
on the chosen fit. 

These error sources were evaluated together by re- 
extracting all the spectra in different manners and perform- 
ing new measures with different fits in the H/j wavelength. 
We repeated 75% of the measurements (326 out of 432), cho- 
sen in a stochastic way in order to avoid selection effects in 
this choice. The distribution of the differences among these 
new measurements and the ones used for this paper give the 
combined extraction and fit error cTgx (Table |3}. We divided 
the spectra in groups with similar S/N, i.e. hot and cool stars 
(Teff ^ 20 000 K) and very hot stars in field B (T^g * 30 000 K). 
In field A, the extraction and fit eiTors for the CCs involving the 
not-summed A14 spectra (much noisier) were evaluated sepa- 
rately and, as expected, they came out much greater than the 
others. 

4.4. Sky line position errors 

The coiTection procedure described in ^3. II introduced an ad- 
ditional error, due to the uncertainty in the measurement of the 
5577 A sky line position. We assumed this error to be equal to 
the dispersion of the corrections around the least-square solu- 
tions when plotted against the slit Y-position, as in Fig.|5] We 
estimated cTsj^y = 0.025 A, as expected for a fit eiTor cr = 0.1 
pixels. A similar value was obtained by analyzing the Gaussian 
fits of the sky lines. Then for the measurements in we have 

cr,ky = 0.025 A = 1.54 km S-' 

to be considered twice (with proper eiTor propagation), as each 
correction involves two independent measurements. 

4.5. Correction to RV measures 

The corrections to the RV measurements described in ^3.2l have 
been an additional source of error. For each pair of frames, we 
obtained the corrections by plotting the differences in stellar 
positions as in Fig. [7] The standard deviations of residuals with 
respect to the least-square solution were calculated for each 
pair, and used as the estimated error o-com introduced by this 



procedure. It changed only slightly from one pair to the other, 
being always 0.7 < crcon < 1.4 km s"'. 

5. Results 

5.1. RV variations with 

We obtained the most accurate and reliable RVs using the Hp 
line. The RV variations from all the possible combinations of 
the observed spectra are listed in Table |3 For each star, Fig.|9] 
shows the absolute value of the measured maximum RV varia- 
tion. 

All the variations are small, usually below 10 km s and 
never further than 3cr from zero; therefore no RV variation can 
be considered significant. The variations are only slightly larger 
in the blue part of the HB, but with larger errors due to decreas- 
ing S/N in the spectra. 

The highest variation plotted in Fig.|5](21.7 km s"') refers 
to star 2823 1 . We consider this datum interesting but particu- 
larly dubious, since it was measured in frame A14 where we 
measured higher variations, due to higher noise of the spectra, 
and star 28231 showed particularly distorted line profiles. This 
variation is not confirmed by the measurements in other wave- 
lengths, and this star shows no other remarkable variation on 
the other nights. Therefore, it is highly probable that the rela- 
tively high RV variation is only due to noise-distorted Hp wing 
profile. 

Table0]also gives the results concerning the frame B 1 3 (un- 
corrected for the star position inside the slit, see ^3.2> . although 
these have been excluded from our analysis and from the plot in 
Fig- El It can be seen that this exclusion represents no substan- 
tial loss of information, since the RVs appear scattered around a 
certain systematic value different from zero, with no larger RV 
variations than the typical 3cr interval. Then even when includ- 
ing the data from the uncorrected frame B13, we can confirm 
our general results. 

5.2. RV variations in other wavelengths 

The results with Hy, when available in the observed spectral 
range, always confirm the ones using Hp. Occasionally the dif- 
ference reaches 10 km s"' (Fig.^J central panel), but this 
agrees with a Gaussian distribution of the differences with a 
dispersion equal to our estimated errors. 

The measurements with weak metallic lines gave very good 
results for cool stars (Fig. ^| upper and bottom panel), be- 
cause of the large number of lines and the good S/N of the 
spectra. The RV variations are always small, and the general 
trend confirms the results we had using Hp (the differences are 
on average below 3 km s"'). For hot stars the results obtained 
using weak lines (mainly He lines) are not very reliable for 
many reasons: they are not uniform, since the number of useful 
lines changes from spectrum to spectrum; the wavelength is not 
clearly determined, so the applied coiTections introduce some 
uncertainty, due to the translations from A to velocity quanti- 
ties; the lines are few and weak, and then easily distorted by 
noise spikes. Still all the RV variations are within the 3<t error 
bars, although the errors are larger than in the previously dis- 
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Fig. 9. Maximum RV variation for each star (in absolute value) measured in wavelength, plotted as a function of Teff. The 
thick errorbars are the Icr errors, while the thinner bars indicate the 3<t. There is no RV variation exceeding the 3cr error bars. 



cussed measurements. For some stars the maximum RV varia- 
tions occasionally reach 30 km s"' (Fig. [TO] bottom panel). It 
is worth noting that these differences always tend toward lim- 
iting high variations in Hj}, never toward emphasizing them, 
and this is also clear from the fact that, in spite of these great 
differences, the maximum RV variations are on absolute value 
that is the same order of magnitude as in or just slightly 
higher. For example, both the stars 28231 and 28947, which 
showed the highest RV variations in H^, show no great vari- 
ation in the weak line measurements. Therefore the measure- 
ments with weak lines confirm the results with H^, and indicate 
that the higher RV variations observed with Hp are simply due 
to random errors. 



5.3. Binary detection probability in our observations 

In order to better understand the significance of our results we 
calculated the probability d of detecting a binary in our obser- 
vations as a function of the period P. In order to relate the period 
and the maximum semiamplitude of the RV curve, we assumed 
a circular orbit and a mass of 0.5 Mq for both components. 
These assumptions are representative of the typical binary sys- 
tems observed in field sdB stars. 

For each value of P we considered 50 possible values of 
vsin(/) (where ; is the inclination of the orbit along the line 
of sight), and 50 possible values of the phase To (equally 
distributed at constant step in the range < Tq < P and 
< V sin(i) < 1). Then, we calculated how many of these 2500 
binary systems would have been detected in our observations, 
defining as "detection" an RV variation of more than 20 km s ' 
(our 3cr for hot stars) between any two of our observed epochs. 
Finally we weighted the probabilities for the three fields by the 



number of hot stars observed in each, in order to derive the 
average detection probability d of our observations. 

The results are shown in Fig. ^2 The probability of detect- 
ing a binary with periods P< 5 days is usually higher than 80%, 
and reaches 90% for periods shorter than one day. There is an 
evident loss of sensitivity around P = 1 day, since this is the 
typical temporal interval between two observing epochs. 

We performed similar calculations also assuming a com- 
panion of 0.1 Mq. This kind of system is a minority among the 
binary population in the field, but t hey do exi st, as for example 
the well-known HW Vir system (Menzies & Marang.1986.) . As 
shown in Fig the detection probability for such systems is 
very low in our survey, with the exception of the shortest (1-2 
days) period binaries. 

5.4. Binary fraction estimate 

If the fraction of binaries in the sample is / the probability of 
detecting Nb binaries out of a sample of N stars is: 



P = 



N! 



(dffHl-df) 



(N-Nb)!Nb! 



N-Nb 



where d is the probability of detection weighted for the period 
distribution of binaries. In our survey Nb=0 and N=18, ther- 
fore, 

P = (1 -df)'l 

Since we found no binaries, the probability P has a unitary 
maximum in f-Q. Then, for increasing values of /, it falls 
rapidly to zero. The exact shape of the function depends on 
the assumed period distribution, which is very uncertain. We 
performed the calculations for two limiting cases, a flat distri- 
bution and a truncated Gaussian centered on P= 1 day and with 
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Fig. 10. Summary of the results in other wavelengths, plotted vs. effective temperature of the stars. Top panel: differences between 
the results in and the weak lines. Central panel: differences between the results in and Hy. Bottom panel: maxima RV 
variations measured with weak lines. The thick error bar indicates the Icr interval, and the thin one is the 3cr. 



log(P/ davs )=l as FWHM, as also assumed by iMaxted et alJ 
(1200 1*) and'N apiwotzki et alJ l EoO ?). This distribution seems to 
follow the observations of field sdB stars, although the known 
data are still too few and this is just a first guess. The results 
on P obtained with the two distributions are extremely simi- 
lar to within 1 %, because the decline of the Gaussian wings in 
the second case compensates for higher (lower) values of the 
sensitivity for short (long) periods, and the weighted mean is 
very similar The value of P reaches 0.05 for /=0.20, in both 
cases. Then we can conclude that the binary fraction among 
EHB stars in NGC6752 is lower then 20% at a confidence level 
of 95% . 

6. Discussion and Conclusions 

The results of Fig.|9lclearly indicate that there is no close bi- 
nary system in our sample. The RV variations observed dur- 
ing the four nights, with different temporal samplings, are all 
within the estimated errors. These errors (a; 4 - 6 km s"') are 



significantly smaller than the radial velocity variations we were 
expecting based on the results for field sdB stars. 

In the compilation bv lMoral~ Rueda et al.1 ( l2003l) . for all 
the 38 field sdBs with known orbital periods, the RV semi- 
amplitude is always greater than 30 km s"'. These binaries 
should have been easily detected in our survey at a 5-6 sigma 
level. If the binary fraction of the EHB st ars in NGC6752 is th e 
same as among the field sdB stars listed in lMaxted et al ] ll200ll) . 
we would expect that 13 ± 2 of the 18 stars with Te >22,000 K 
in our sam ple should be close bi naries; given the period dis- 
tribution of iMorales-Rueda et all 12003) . 80% of them should 
have a period P < 5 days and, therefore, be easily detected with 
our temporal sampling. Indeed, as shown in the previous sec- 
tion, our survey should be able to detect 85% of close binaries 
with P < 5 days of the kind found among the field sdB stars. 
I.e., if the binary fraction among the hot cluster EHB stars were 
the same as among field sdB stars, we would have expected to 
find 9 binaries with RV variations > 20 km s ' . 
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Fig. 11. Probability of detecting a binary system in our obser- 
vations as a function of its period P. The solid line indicates the 
probability calculated for an sdB star and a companion of 0.5 
M0, the dotted line for a companion of 0.1 M0, assuming an 
EHB star of 0.5 M0 in both cases. 



We found no significant RV variation at 3 sigma level (15- 
18 km s^'). Only 18 measurements showed a variation exceed- 
ing 2 sigmas (10-12 km s"'). However, assuming a Gaussian 
distribution for our errors, we would expect that 19.6 mea- 
sures among our 432 RV difference estimates should exceed 
the 2 sigma level; therefore, these variations too cannot be con- 
sidered significant. As demonstrated in Sect. 15.41 these results 
imply that the fraction of close binaries among EHB stars in 
NGC6752 is < 20%, at a 95% confidence leve l This result 
is in sharp contrast with that of IPeterson et al ] ll2002h . who 
concluded, from the measurement of radial velocity variations 
of 30 EHB sta rs in NGC6752 , that the majority of them are 
close binaries. IPeterson et al ] iE002.^ unfortunately have not 
published more details about their investigation so far, and 
therefore, it remains to be established whether the d i fiferent 
conclusions are due to the lower S/N o f 'Peter son etaP (l2002h 
data (who used the multifiber Hydra spectrograph mounted on 
a 4m size telescope) or to some physical reasons. Interestingly 
enough, apparently, our sample of stars is different from their 
sample. In particular, according to Peterson (priv. comm.), their 
EHB stars are on average located in a less crowded, more ex- 
ternal environment than ours (though our stars are evenly dis- 
tributed from ~ 100" to ~ 500" from the cluster center, i.e. out 
to about 3.5 half-mass radii, see Fig.|2j. We have only five stars 
in c ommon, of wh i ch on ly two are EHB stars and for which 
also lPeterson et al ] |l2002ll do not find evidence of radial veloc- 
ity variations (Peterson, priv. comm.). 

In summary, we are forced to conclude that most (> 80%) 
of the hot EHB stars in NGC6752 are not close binaries. This 



result poses a number of problems. First of all, it implies that 
the mechanism for the formation of EHB stars in GCs does not 
involve a large envelope mass loss enhanced by tidal interac- 
tion within a close binary at variance with what happens for 
the field sdBs. Or, at least, this cannot be the dominant forma- 
tion mechanism. In a sense, this result is not completely un- 
expected. The typical companions of field sdB s tars are r ela- 
tively massive stars (0.3-0.5 Mq, see iMaxted et aDl2004l) . If 
the same binary systems are the source of the EHB stars in 
GCs, they should originally have had a mass that is on aver- 
age larger (> IM0) than the mass of a single star at the TO 
(~ 0.8 M0) and, therefore, should be more concentrated than 
the upper MS/RGB stars because of mass segregation. Despite 
a number of attempts, there is no compelling evidence at the 
moment that EHB stars are, on average, more concentrated than 
cooler HB stars in the same clusters, as we would expect if they 
were members of a binary system with similar properties to the 
binaries producing field sdBs, and as we find, e.g., for the blue 
stragglers in all the obse rved GCs (Piotto et al. 2004). Instead, 
recent surveys (e.g. by Bedin et al. 2000, in NGC2808, a pro- 
totype EHB cluster) have shown that no clear radial gradient of 
the EHB stars can be detected. This empirical evidence, cou- 
pled with the results of the present paper, implies that: either 1) 
the companions of the EHB stars have very low masses (so that 
the sum of the main sequence mass of the EHB stars and that 
of the companion do not differ much from the mass of a typical 
MS star), or 2) cluster EHB stars all have periods much longer 
than 5 days, or 3) cluster EHB stars are not binaries. 

In any case, we are forced to conclude that in GCs there are 
diff'erent formation mechanisms of EHB stars with respect to 
the main formation channels of field sdBs. It is possible that the 
dynamical evolution of GCs have removed the primordial bina- 
ries able to produce sdB stars, at least in the inner part o f the 
cluster. This is an interesting possibility. Indeed, Piotto e t alJ 
(12004ft have shown that the frequency of blue stragglers (BSs) 
in the cores of GCs is significantly smaller than the frequency 
of field BSs, and that there is an anticor relation between the 
frequency of BSs and the GC total mass. IPavies et al.l (l2004l) 
suggest that this is likely due to the fact that in more massive 
clusters, where the encounter probability is higher, binary evo- 
lut ion is accelerated by si milar mechanisms to those proposed 
bv lPavies & Hansenl ( ll99 8) to explain the formation of pulsars 
in GCs. In this way, the frequency of BSs (from primordial bi- 
naries) could have been much higher in the past, but now there 
are much fewer binaries to produce them via the merger chan- 
nel. However, the GC environment can favor different forma- 
tion mechanisms fo r the EHB stars. Ag ain, something similar 
happens for the BSs. lPiotto et alJ ( [20041) show that the luminos- 
ity functions of the BSs in the core of massive clusters diff'er 
significantly from those in less massive ones, and the former is 
compatible with BSs mainly formed via colUsion of stars (more 
probable in more massiv e clusters). 

Interestingly enough, Rec io-Blanco et aP ( l2005l) show that 
more massive clusters also show more extended HBs. If the 
IPeterson et al.l J2002) results are confirmed by a more extended 
survey, i.e. if at least a fraction of the EHB stars in the out- 
skirts of the clusters are indeed close binaries (at variance with 
what we find in the more internal regions), we will have fur- 
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ther evidence of the effect of the environment on the formation 
of EHB (or sdB) stars. This would also explain the sharp dif- 
ference between the frequency of close binaries among clus- 
ter EHB st ars and field sdBs. T he results of the present pa- 
per and of IPeterson et alJ (l2002l) might suggest that two dif- 
ferent channels of production of EHB stars are at work, at 
least in some GCs. Again, it is tempting to recall that there 
is compelling evidence of two channels of BS production in 
the same GC. It has been shown that the radial distribution of 
BSs in M3 ( Ferraro et al. 1997) and in 47Tuc ( Ferraro et al. 
I2OO4I) are clearly bimodal and that the BS properties in the 
inner regions and in the outskirts of this cluster are also dif- 
ferent. [Mapelli et al. (2004) show that the radial distribution of 
BS in 47Tuc supports the hypothesis that the BSs produced in 
the center are caused mainly by collisions, while the BS in the 
cluster outskirts (for r > 30 core radii, i.e. for distances from 
the center greater than 4 half mass radii) must all result from 
the merge of primordial binaries, which survived because of 
the much slower dynamical evolution in the external part of the 
cluster A similar scenario could also be at work for EHB stars, 
with EHB stars in the external part of the cluster coming from 
the close binary-evolution channel suggested for the production 
of field sdB stars. It would be interesting to use the appropriate 
dynamical and stellar evolution models to investigate the pos- 
sibility that EHB stars in GCs can be formed preferentially by 
stellar collisi ons in their mo r e cent ral parts. Surely, our results 
and those of IPeterson et al ] (l2002h call for a much more ex- 
tended observing campaign to search for close binaries among 
EHB stars in this and other GCs. 

Another interesting possibility that needs to be (observa- 
tionally) explored is that cluster EHB stars are produced in 
binaries with significantly longer periods or with significantly 
less massive companion than in field sdBs. As shown in Figll II 
our survey could find only binaries with relatively large-mass 
companions, and as it only has a 4-day temporal coverage, it is 
totally insensitive to t hose binaries with periods longer than 10 
days. We note that the lHan et al ] ll2003l) models for the produc- 
tion of sdB stars predict a large number of long period (~ 100 
days) binaries. A follow-up observing campaign is needed to 
verify if there are binaries like these among cluster EHB stars; 
in any case, long-period binaries seem to be a minority among 
field sdBs, so that the question to why field and cluster EHB 
stars should form through different mechanisms remains unan- 
swered. 

Of course, it is equally possible that dynamical evolution 
and binaries have nothing to do with the formation of EHB 
stars in GCs. Unlike field sdBs, which can have progenitors 
with rather different masses, all the cluster EHB stars start their 
evolution from a well defined TO mass (~ 0.8 M0, slightly de- 
pending on cluster age and metallicity), because of the very 
small (compatible with zero) age and metallicity dispersion 
among stars in a given GC. There must be an extremely well- 
tuned (unlikely?) mass-loss mechanism to produce a star with 
the very small core and envelope mass dispersion of the clus- 
ter EHB stars. Other explanations are possible. As discussed 
in the introduction, recent results on very mass ive clusters like 
u) Cen tauri (Piotto et al..2005.) and NGC2808 (iP' Antona et alJ 
l2005h seem to indicate that the EHBs in GCs could represent 



the evolved population of a second generation of stars formed 
by material enriched in He because of the pollution by the 
ejecta of SNe and/or intermediate mass AGBs from a first gen- 
eration of stars. In this hypothesis, it is possible that more mas- 
sive clusters are better able to keep part of the ejecta, ex plain- 
ing the correlation found by iRecio-Blanco et al ] (I2OO5I) . One 
should keep in mind, however, that not only the total mass, but 
also the concentration of a globular cluster determine the depth 
of its gravitational potential and thereby its ability to retain en- 
riched material. 

There is an interesting point to add before conclud- 
ing this discussion. As noted in the introductory section, 
Napiwotzki et al. (2004) found a much lower (42%) binary fre- 
quency among field sdBs than did lMaxted et al.l ( l200ll) . On one 
hand, even with this low binary frequency, we would expect to 
find 5 binaries arnong ou r sampled 18 EHB stars in NGC6752. 
iNapiwotzki et al.l ll2004l) suggest that the difference in mean 
apparent magnitude between their sample and the one of 
Maxted et al. (2001) may imply a difference in the populations 
sampled by the two surveys. Although this has not yet been 
proven, attributing the observed difference in binary fraction to 
a population difference would mean that binaries are substan- 
tially less frequent among thick disk-halo sdBs, and it would 
imply a trend with age and/or metallicity. On average, our clus- 
ter EHB stars are expected to be older, and possibly more metal 
poor than the Napiwotzki et al. sample, and, apparently, the 
fraction of binaries among them is even smaller than among the 
Napiwotzki et al. sample. It is currently unclear if these differ- 
ences are related to the different environments, ages, or metal- 
licities, or to a combination of these parameters. Supporting 
evidence for the possible influence of abundance differences 
comes from the obs erved abundance a nomalies in GC red gi- 
ants ( Catelani2005tiGratton et alJ20 04'. and references therein), 
which are not seen in field red giants (Gratton et al. 2000 ). As 
suggested by niany in the Uterature (e g.lVandenberg & SmithI 
'l988l lSweigartlll997tlD^ntona etalJl2002h . these abundance 
anomalies (either primordial, or due to mixing effects) can af- 
fect the RGB evolution and, consequently, have an impact on 
the properties of HB stars including temperatures, luminosities, 
gravities, and pulsation characteristics. 
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Table 4. RV variations for target stars. 
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0.4 


-3.8 




3.8 ±4.7 


4.8 ±4.5 




14.3 - 


-3.2 ±4.3 


-1.8 ±4.2 




4.2 


-4.7 ± 3.9 


9b 


-2.7 ± 2.5 


-3.3 


-6.3 




11.0 ±3.9 


-8.6 ±3.7 




4.0 - 


-8.4 ±3.6 


-5.8 ±3.7 




2.3 


2.2 ± 3.4 


10b 


-2.1 ± 3.2 


2.9 


0.2 




-1.5 ±4.6 


1.5 ±4.5 




6.8 - 


-3.8 ±4.2 


-0.9 ± 4.4 




2.1 


-1.7 ±4.6 


lib 


-1.2 ± 2.4 


1.3 


-3.0 




1.5 ± 3.9 


2.8 ±4.0 




11.6 


2.9 ±3.6 


4.0 ±3.6 




9.3 


1.4 ±3.8 


12b 


-1.1 ± 3.5 


0.7 


-4.7 




-0.5 ± 4.6 


-0.4 ± 4.5 




9.7 


1.1 ±4.9 


3.9 ±4.9 




10.9 


2.5 ± 4.7 


13b 


1.1 ± 3.8 


3.6 


-2.6 




3.2 + 5.1 


1.3 ±5.3 




9.4 - 


-3.8 + 4.5 


-5.4 ± 4.5 




-0.4 


-7.4 ± 5.4 


14b 


0.9 ± 3.2 


6.8 


0.9 




-1.0 + 4.3 


-2.2 ± 4.5 




3.3 


3.0 ±4.2 


2.1+4.1 




3.6 


2.8 + 4.1 






star 




14a vs 13 


14b vs 13 14b vs 14a 15 vs 13 


15 vs 14a 15 vs 14b 












Ic 




-0.3 ± 3.7 


-2.3 


±3.8 -2.0 


± 3.7 -0.9 


±3.8 


-0.6 ± 


4.0 1.4 


±3.4 












2c 




0.1 ±3.7 


-1.9 


±4.0 -1.7 


± 3.7 -4.3 


±3.7 


-4.2 ± 


3.7 -2.2 


± 3.5 












3c 




-1.3 ±3.8 


-2.0 


± 3.9 -0.4 


±3.6 -2.1 


±3.8 


-0.8 ± 


3.8 -0.1 


±3.5 












4c 




1.4 + 3.8 - 


-0.2 


+ 4.0 -1.5 


+ 3.6 -0.1 


+ 3.8 


-1.6 ± 


3.9 0.0 


±3.4 












5c 




2.5 ± 3.7 


4.8 


±3.7 1.7 


±3.6 3.2 


±4.0 


0.3 ± 


3.7 -2.0 


± 3.4 












6c 




1.4 ±4.0 


0.4 


±4.2 -0.3 


± 3.6 0.3 


±3.8 


-0.9 ± 


3.7 -0.6 


± 3.4 












7c 




3.4 ±3.9 - 


-1.4 


± 3.8 -4.9 


± 3.6 -4.3 


±3.8 


-7.5 ± 


3.7 -2.8 


± 3.7 












8c 




-1.9 ±4.5 


3.9 


±4.9 5.9 


± 4.5 6.2 


±4.5 


6.8 ± 


4.5 -0.8 


±4.5 












9c 




5.4 ±3.8 


1.8 


±3.7 -3.6 


±3.7 3.2 


±3.7 


-2.1 + 


3.7 1.4 


±3.5 












10c 




4.2 ±4.7 


0.5 


±5.0 -2.3 


±4.5 -3.9 


±4.5 


-5.5 ± 


4.5 -4.5 


±4.2 












11c 




-0.7 ± 3.8 


0.4 


±3.7 1.1 


±3.6 3.1 


±3.8 


3.8 ± 


3.7 2.9 


±3.3 












12c 




-1.4 ±4.5 - 


-7.0 


± 4.5 -4.4 


±4.5 -4.0 


±4.7 


-0.8 ± 


4.5 3.2 


±4.5 












13c 




-1.4 ±3.8 


0.4 


±3.7 1.9 


± 3.7 0.6 


±3.7 


1.6 ± 


3.7 0.3 


±3.3 












15c 




2.7 ±3.9 


2.0 


±3.7 -0.8 


±3.7 1.1 


±3.7 


-1.6 ± 


4.0 -0.9 


±3.4 












16c 


10.7 ± 4.4 


3.5 


±4.6 -5.1 


± 4.5 5.0 


±4.5 


-3.9 ± 


4.4 3.2 


±4.2 












17c 




-0.2 ±3.7 


0.5 


± 3.7 0.6 


±3.6 -1.3 


±3.8 


-0.9 ± 


3.7 -1.6 


±3.3 












18c 




2.6 + 3.7 


4.1 


+ 3.7 1.5 


+ 3.6 1.2 


+ 3.8 


-1.9 + 


3.8 -3.4 


±3.4 









